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Introduction 

Several  modem  space  propulsion  devices  exploit  the  acceleration  of  propellant  via  the 
production  of  a  magnetized  plasma,  some,  in  conditions  where  there  are  crossed  electric  and 
magnetic  fields  (e.g.,  Hall  thrusters,  Cusped  Field  Thrusters).  These  thrusters  share  several  features 
that  are  in  common  with  other  magnetized  cross-field  plasma  sources  such  as  magnetron  and 
Penning  discharges.  Magnetized  plasmas  in  cross-field  configurations  exhibit  complex  nonlinear 
behavior  resulting  in  variety  of  turbulent  and  coherent  fluctuations  that  can  critically  affect  thruster 
operation  and  performance.  Rotating  structures  have  been  observed  in  a  variety  of  laboratory  and 
technological  plasma  devices  with  magnetized  electrons  and  non-magnetized  or  weakly 
magnetized  ions.  These  oscillations  usually  have  low  mode  number  with  a  characteristic  frequency 
that  lies  between  the  electron  and  ion  gyro-frequencies,  i.e.,  Qce «  ©  <  Qd .  Experimental  studies 
of  Hall  thrusters  and  Penning  discharges  demonstrated  that  low  frequency  (10’s  kHz)  spoke 
oscillations  contribute  to  anomalous  electron  cross-field  transport.  Studies  carried  out  in  smaller 
devices,  such  as  mesoscale  magnetrons,  place  these  coherent  fluctuations  at  much  higher 
frequencies  (100’s  of  kHz),  and  also  contribute  substantially  to  the  electron  transport.  For  Hall 
thrusters,  electron  cross-field  transport  is  of  practical  importance  because  it  diminishes  the  thrust 
efficiency.  Fluctuation-induced  transport  reduces  the  local  electric  field,  potentially  leading  to 
increased  plume  divergence  and  thruster  wall  erosion.  This  program  is  to  better  understand  the 
physics  behind  these  fluctuations  and  their  relevance  to  propulsion  devices.  Our  experiments  are 
primarily  carried  out  on  Penning  discharges  (Princeton  Plasma  Physics  Laboratory  [PPPL]  and 
Princeton  University),  and  on  magnetron  discharges  (Stanford),  with  some  studies  for  comparisons 
to  cylindrical  Hall  thrusters  (PPPL  and  Princeton  University)  and  cylindrical  cusped  field  thrusters 
(Stanford),  referred  to  here  as  CHT  and  CCFT  for  short.  Both  institutions  have  also  been  actively 
developing  theory  and  simulations  of  the  development  of  these  coherent  structures. 
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Research  Background 

Magnetized  plasmas  in  cross-field  discharge  devices  exhibit  complex  nonlinear  behavior 
resulting  in  a  variety  of  turbulent  fluctuations  and  structures  that  critically  affect  operation  and 
performance  of  these  devices.  ExB  rotating  structures  sometimes  referred  to  as  rotating  spoke 
oscillations  have  been  observed  in  a  variety  of  laboratory  and  technologically  important  plasma 
devices  with  magnetized  electrons  and  non-magnetized  or  weakly  magnetized  ions.  These  devices 
include  Hall  thrusters,  Penning  discharges,  and  sputtering  magnetrons.  These  oscillations  usually 
have  low  mode  number  with  a  characteristic  frequency  of  Qce  «  co  <  Qd,  where  Qce  and  Qd  are 
the  gyrofrequency  of  electrons  and  ions,  respectively.  Recent  experimental  studies  of  Hall  thrusters 
and  Penning  discharges  demonstrated  that  low  frequency  (1-10’s  of  kHz)  spoke  oscillations  are 
responsible  for  anomalous  electron  cross-field  transport.  In  the  magnetron  thruster  studies  with 
segmented  anodes  (yet  unpublished,  see  below),  these  instabilities,  which  appear  to  be  gradient- 
drift  driven  in  this  device,  are  also  responsible  for  the  cross-field  discharge  current).  For  Hall 
thrusters,  electron  cross-field  transport  is  of  practical  importance  because  it  diminishes  the  thruster 
efficiency,  which  is  the  ratio  of  the  plasma  jet  power  to  the  input  electric  power.  If  the  electric  field 
decreases  with  increased  cross-field  transport,  it  can  cause  plume  divergence  and  thruster  wall 
erosion.  Our  recent  measurements  in  Cylindrical  Hall  thrusters  (CHT)  (Fig.  la)  and  Penning 
discharges  (Fig.  lb)  revealed  a  strong  dependence  of  spoke  oscillations  on  the  input  discharge 
parameters,  including  gas,  gas  pressure,  the  magnetic  field,  and  the  electron  injection  from  the 
cathode[l,  2],  Similar  observations  have  been  reported  for  conventional  state-of-the-art  Hall 
thrusters. 

The  moving  spoke  is  a  plasma  non-uniformity  that  rotates  usually  in  the  ExB  direction.  The 
spoke  can  appear  in  different  modes,  but  for  CHTs,  Penning  discharges  and  cylindrical 
magnetrons,  the  most  common  is  the  m  =1  mode  [1-3]  (Fig.  2).  In  CHT  and  Penning  discharges, 
the  spoke  was  observed  with  a  high-speed  CCD  camera  and  Langmuir  probe  measurements  [3], 

a)  b) 
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Figure  1.  Schematics  of  ExB  configurations  of  the  cylindrical  Hall  thruster  (a)  and  the  Penning  plasma 
discharge  (b)  [1,2]. 


The  spoke  appeared  on  the  camera  as  a  spot  of  increased  light  emission  that  rotates  in  the 
azimuthal  direction  with  the  translation  speed  of  an  order  of  magnitude  slower  than  the  local  ExB 
drift  speed.  Similar  observations  were  also  reported  for  the  annular  Hall  thrusters  [4]  and  other 
plasma  devices.  On  probes  immersed  in  CHT  and  Penning  plasmas,  the  spoke  manifests  itself 
through  azimuthal  oscillations  of  plasma  potential,  electron  temperature  and  plasma  density. 

Naturally,  the  spoke  is  a  consequence  of  one  or  coupled  multiple  instabilities  of  the  cross-field 
plasma  discharge.  Our  recent  studies  on  Penning  discharges  provided  new  insights  on  the 
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fundamental  role  of  gradient  drift  (modified  Simon-Hon  type)  instabilities  driven  by  gradients  of 
plasma  density,  temperature  and  magnetic  field  [5],  Ionization  instabilities  may  be  coupled  to  this 
electrostatic  instability  giving  raise  the  formation  of  rotating  plasma  structures.  We  have  identified 
in  the  experiment  and  simulations  a  strong  connection  between  these  plasma  structures  and 
electron  transport.  This  connection  is  responsible  for  the  mode  transition,  facility  effects,  and 
effects  of  the  magnetic  field,  effects  of  device  (e.g.  thruster)  geometry  and  materials  on  plasma 
structures. 


a) 


2.6  cm 

1 - 1 

* 

< 

Time:  (Vs 

Time:  17»a 

Time:  34ns 

¥ 

* 

Time:  51|is 

Time:  68|is 

Time:  85fis 

Figure  2.  ExB  rotating  spoke  in  CHT  thruster  (a)  and  Penning  discharge  (b). 


The  overall  physics  which  we  have  discerned  from  our  integrated  experimental,  simulations 
and  theoretical  efforts  (described  below)  is  as  follows:  small-scale  phenomena  such  as 
microturbulence  is  responsible  for  anomalous  transport,  which  determines  both  the  electric  field 
and  pressure  gradients,  which  in  turn,  can  induce  large-scale  gradient-driven  instabilities  forming 
plasma  strictures  such  as  spokes.  The  dynamics  of  the  plasma  structure  is  also  determined  by  a 
relatively  slow  rotation  of  ions  due  to  a  partial  magnetization  or  momentum  conservation,  and 
charge  exchange  collisions. 

Understating  of  micro  turbulence  at  the  scale  comparable  to  the  electron  gyro  radius  and  its 
effects,  including  the  formation  of  the  plasma  structures  requires  a  kinetic  description  supported 
by  theory  and  validated  by  measurements  of  kinetic  properties,  including  energy  distribution 
functions  of  electrons  and  ions.  We  believe  that  it  is  important  to  conduct  such  simulations  and 
measurements  with  a  focus  on  three-dimensional  ExB  plasma  configurations  with  non-uniform 
magnetic  fields,  operating  under  the  realistic  plasma  conditions  that  are  encountered  in  plasma 
propulsion  devices  and  anticipate  this  as  a  major  component  of  any  continued  research  project. 

Our  studies  have  taught  us  that  the  electron  plasma  response  is  crucially  sensitive  to  the 
dynamics  along  the  magnetic  field  direction  due  to  high  velocity  of  electron  thermal  motion.  On 
the  contrary,  the  motion  across  the  magnetic  field  occurs  on  a  slower  time  scale  due  to  ExB  drift 
and  drifts  due  to  magnetic  field  inhomogeneity  and  even  slower  drifts  due  to  electron  inertia  or 
collisions.  The  ratio  of  the  transverse  to  the  parallel  times  scale  can  change  the  conditions/nature 
of  the  instabilities  in  a  significant  manner  and  affects  anisotropy  of  electron  velocity  distribution 
function.  Such  changes  can  be  induced  by  boundary  conditions,  by  changing  the  geometry  of  the 
plasma  set  up,  by  driving  electron  fluxes  along  the  magnetic  field  between  electron  emitting  and 
electron  absorbing  walls  applied  to  control  density  gradients  across  the  magnetic  field,  or  by 
shortening  the  plasma  along  the  magnetic  field  and  thereby  affecting  gradient-drift  instabilities,  or 
by  exciting  ion-acoustic  instability  along  the  magnetic  field,  etc..  The  first  approach  resembles  so- 
called  segmented  Hall  thruster  concepts  proposed  at  PPPL  more  than  a  decade  ago  (Fig.  3)  [6,7], 
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but  never  considered  or  studied  in  terms  of  its  effect  on  plasma  instabilities  and  plasma  structures. 
The  results  of  our  recent  study  indicated  that  changes  in  boundary  can  strongly  affect  anomalous 
transport  [8]. 


Accomplishments 


Experiments  on  Penning  Discharges 

Experiments  at  PPPL  were  conducted  primarily  using  the  Penning  discharge  system  (Fig.  lb). 
With  a  flexible  and  versatile  design,  this  system  has  a  convenient  access  for  various  diagnostics 
and  can  operate  with  well-controlled  plasma  parameters.  In  its  operation,  a  dc  voltage  of  20-200 
V  is  applied  between  the  RF  plasma  cathode  and  the  anode-chamber.  The  discharge  can  be 
operated  with  either  a  uniform  or  non-uniform  magnetic  field  by  varying  the  directions  of  the 
currents  in  the  electromagnet  coils  or  by  powering  coils  with  unequal  currents.  In  all  previous 
experiments,  the  Penning  system  was  operated  in  a  Helmholtz-like  configuration  with  co-directed 
equal  currents  in  both  coils  to  produce  a  relatively  uniform  magnetic  field.  When  operating  over  a 
relatively  broad  range  of  magnetic  fields,  of  30-500  Gauss,  the  Penning  discharge  sustains  an 
efficient  ionization  of  xenon  with  plasma  properties  (Fig.  4)  [2]  comparable  to  the  plasma 
properties  in  the  near  anode  region  of  the  Hall  thrusters  where  strong  rotating  spoke  oscillations 
are  usually  observed  [9]. 
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Figure  3.  A  PPPL  concept  of  segmented  electrode 
Hall  thruster  [7],  which  can  be  adapted  for 
controlling  of  gradient  driven  instabilities  in  the 
proposed  research. 


Our  measurements  have  revealed  the 
presence  of  a  spoke  under  almost  all  conditions 
of  the  magnetic  field  strength.  These  spoke 
oscillations  are  found  to  depend  strongly  on  the 
input  discharge  parameters,  including  gas 
pressure  and  applied  magnetic  field  (Fig.  4d). 
The  gas  type  (Xe,  Ar)  was  also  shown  to  affect 
the  spoke  oscillation  (Fig.  5).  These 
measurements  revealed  that  the  spoke 
frequency  follows  the  scaling/oc  \/px(B/m)m, 
where  B  is  the  applied  magnetic  field,  m  is  the 
ion  mass,  and  p  is  the  pressure  [10],  An 
increase  in  the  gas  pressure  above  10'4  to  10'3 
Torr  leads  to  the  suppression  of  the  spoke 
oscillations.  This  scaling  is  different  from  the  scaling  observed  in  the  cylindrical  Hall  thrusters 
(where /oc  (p/B)2  x  (l/m)m)  [11]  and  magnetron  discharge  (where/  oc  p/B2  [10]).  These  diverse 
experimental  scaling  relationships  describing  the  spoke  dynamics  in  different  plasma  devices 
present  a  challenge  to  theories  of  rotating  plasma  structures.  Previous  studies  pointed  to  a  Modified 
Simon-Hoh  Instability  (MSHI)  as  one  of  the  instabilities  which  could  be  responsible  for  the  spoke 
formation  in  Hall  thrusters.  This  electrostatic  instability  was  proposed  by  Sakawa  et  al  [12]  for  a 
Penning-type  ExB  system  with  density  gradients  and  the  electric  field  perpendicular  to  the  uniform 
magnetic  field.  We  showed  that  the  MSHI  occurs  when  E0/B  >  /2CiLn,  where  E0  is  the  applied 
electric  field,  Ln  is  the  characteristic  scale  of  the  density  gradient,  and  E2ci  is  ion  cyclotron 
frequency.  According  to  Sakawa,  the  rotation  frequency  of  MSHI-induced  plasma  structures  is/ 
oc  (E/m)1/2.  This  theoretical  result  implies  that  the  spoke  should  not  be  affected  by  the  magnetic 
field  or  the  background  gas  pressure  if  the  radial  electric  field  is  constant.  This  contradicts  our 
experimental  findings.  However,  our  experiments  have  revealed  that  an  increase  in  the  gas  pressure 
decreases  the  value  of  the  electric  field  (Fig.  4a)  whereas  an  increase  in  the  magnetic  field  leads  to 
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an  increase  in  the  electric  field.  Thus,  changes  in  these  discharge  parameters  cause  changes  in  the 
electric  field.  According  to  MHSI  theory  and  scaling,  changes  in  the  electric  field  could  affect  this 
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Figure  4.  Gas  pressure  effects  on  the  magnetized  plasma  in  the  Penning  discharge.  Results  of  electrostatic 
probe  measurements  of  plasma  properties  across  the  magnetic  field  in  the  Penning  discharge:  plasma 
potential  (a),  rotating  spoke  frequency  obtained  from  high  speed  imaging  for  different  magnetic  fields. 

instability.  The  suppression  of  the  spoke  with  increased  pressure  correlates  with  a  strong  reduction 
of  the  5-field  (to  nearly  zero).  Changes  in  the  spoke  frequency  with  5-field  correlate  with  changes 
of  E/B  in  combination. 

Based  on  results  of  our  PIC  simulations  and  experimental  data  (discussed  below),  we  find  that 
to  explain  the  spoke’s  behavior  it  is  necessary  to  self-consistently  predict  the  electric  field  as  a 
function  of  discharge  parameters.  This  is  a  rather  complicated  task  because  the  electric  field  itself 
depends  on  the  anomalous  transport,  which  depends  on  the  electric  field.  This  makes  the  problem 
strongly  nonlinear.  Anomalous  transport  is  determined  not  only  by  the  spoke  but  also  by  small- 
scale  fluctuations  that  may  have  different  dependencies  on  the  magnetic  field  and  pressure  and  that 
may  vary  between  ExB  devices.  In  fact,  an  analysis  of  our  experimental  results  using  Ohm’s  law 

across  the  magnetic  field  j  =  -^r  (eE  —  Te  ^n/ngr^  suggests  that  the  electron  transport  is 

anomalously  high  with  veff  /vcoi  ~  10-100  when  compared  to  classical  collisionally-driven 
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Figure  5.  Local  power  spectrum  of  the  normalized  light  intensity  from  fast  camera  images  of  the  spoke 
in  the  Pennine  discharge  for  a  given  magnetic  field  of  50  Gauss  and  different  gases,  argon  and  xenon. 
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transport.  Here,  veff  is  the  effective  electron  collision  frequency  and  vcoi  is  the  collision  frequency 
of  electrons  with  xenon  neutral  atoms  and  ions.  Remarkably,  we  find  these  results  to  be  in  fair 
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Figure  6.  Effective  Hall  parameter  a>ec/ veff  from  (a)  probe  measurements  and  (b)  PIC  simulations  at  the 
B-field  of  50  Gauss  and  two  different  xenon  gas  pressures.  Table  shows  effect  of  ion  mass  on  anomalous 
transport  and  the  effective  Hall  parameter. 

agreement  with  PIC  simulations  (Fig.  6). 

Theory  and  Simulations  Related  to  Penning  Discharge  Studies 

The  commercial  Particle-In-Cell  (PIC)  code  LSP  has  been  extensively  modified,  benchmarked 
and  validated  for  low-temperature-plasma  (LTP)  applications  at  PPPL.  The  standard  version  of 
LSP  can  perform  electrostatic  (ES)  and  electromagnetic  (EM)  simulations  in  one,  two  or  three 
dimensions  in  Cartesian,  cylindrical  or  spherical  coordinates.  It  is  a  versatile  and  capable  code,  but 
has  not  been  extensively  used  for  LTP  applications.  While  we  perform  simulations  some 
deficiencies  were  found  and  remedied.  In  two-  and  three-dimensional  simulations,  we  were 
sometimes  unable  to  get  the  existing  iterative  ES  field  solver  to  converge  the  residual  error  to  a 
sufficiently  small  value  to  avoid  unphysical  effects.  We  therefore  took  advantage  of  LSP  flexibility 
to  implement  different  ES  field  solver.  The  LSP  PETSc  interface  was  updated  to  support  the  latest 
version  of  Argonne’s  PETSc  solver  library  [13].  We  chose  two  direct  ES  field  solvers  that  were 
implemented  and  have  been  exclusively  used  for  our  two-  and  three-dimensional  simulations  in 
both  Cartesian  and  cylindrical  coordinates.  The  direct  ES  solver  based  on  the  SuperLU  library 
(called  through  PETSc)  [Error!  Bookmark  not  defined.]  is  very  robust  and  memory  efficient, 
whereas  the  one  based  on  the  MUMPS  library  is  faster  and  more  scalable  [Error!  Bookmark  not 
defined.] .  The  existing  external-circuit  model  was  found  to  be  numerically  unstable  due  to 
neglecting  the  sheath  capacitance  and  a  new  algorithm  was  developed  and  implemented.  We  also 
had  to  fix  injection  models  in  LSP.  The  model  for  direct  injection  of  a  particle  in  the  simulation 
domain  as  well  as  model  for  secondary  emission  used  fractional-weight  particles.  We  found  these 
models  both  numerically  inefficient  and  difficult  to  use  due  to  large  production  of  particles  with 
fractional  weights.  We  therefore  implemented  injection  models  that  create  particles  using  a 
rejection  method  [14].  With  these  new  injection  models,  all  particle  of  a  species  can  be  guaranteed 
to  have  the  same  numerical  weight,  regardless  of  when  and  how  they  were  created.  Most  of  the 
collision  models  were  also  rewritten  to  use  state-of-the-art  algorithms,  including  charge  exchange 
and  anisotropy  models  for  elastic  and  inelastic  collisions,  as  well  as  for  ionization  [15].  The 
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modifications  to  the  standard  LSP  code  were  deemed  significant  enough  to  justify  using  the  name 
PPPL-LSP  for  the  improved  version  of  the  code. 

Benchmarking  and  Validation 

The  one-dimensional  ES  PIC  code  EDIPIC  was  used  to  benchmark  LSP  for  a  glow -discharge 
validation  for  both  codes,  primarily  their  collision  models  by  comparison  of  the  electric-field 
profile  for  the  two  codes,  as  well  as  for  the  experiment  [16].  To  get  good  agreement  with  a  glow 
discharge,  we  found  that  a  PIC  code  has  to  have  very  accurate  collision  models.  A  paper  on  this 
important  benchmarking/validation  effort  is  under  preparation  [17].  Further  validation  of  PPPL- 
LSP  was  done  for  a  diocotron  instability  frequently  observed  in  Malmberg  traps  or  beam 
experiments  [18].  This  proves  that  implemented  ES  solver  can  resolve  fine  structures. 

Simulations  of  Penning  discharge 

Simulations  were  set  up  with  PPPL-LSP  code  to  investigate  coherent  structures  and  anomalous 
transport  in  a  Penning  discharge.  The  PIC  simulations  resolve  thin  Debye  radius  and  fast  electron 
plasma  and  cyclotron  frequencies.  Therefore  to  reduce  the  simulation  time,  the  device  size  was 
reduced  by  and  order  of  magnitude.  We  made  sure  that  electrons  are  magnetized  and  ions  are  not 
(their  corresponding  gyroradius  smaller  and  larger  than  discharge  radius,  respectively).  To  avoid 
numerical  problems  on  the  axis,  Cartesian  coordinates  were  used  instead  of  cylindrical  [19].  A 
homogeneous  magnetic  field  was  applied  with  a  typical  value  of  100  Gauss.  Similarly  to 
experimental  set  up,  a  narrow  electron  beam  was  injected  along  the  axis.  To  delineate  ionization 
and  transport  mechanisms  of  coherent  structures  formations,  injected  beam  energy  was  below  the 
ionization  threshold  in  initial  simulations.  The  plasma  was  surrounded  by  a  grounded  wall  (anode). 
For  most  of  the  simulations,  the  simulated  device  was  filled  with  helium  gas  at  0.2  mTorr.  In  this 
configuration,  the  electric  field  is  not  externally  controlled  as  in  magnetron  discharges  but  is 
established  self-consistently  by  necessity  to  conduct  radial  current  flows  radially  across  the 
magnetic  field,  which  is  needed  to  balance  the  injected  axial  current.  This  configuration 
corresponds  to  experimental  set  up  where  voltage  is  established  as  a  function  of  the  cathode 
current,  and  the  radial  current  is  maintained  either  by  an  electric  field  or  by  a  pressure  gradient 
with  the  electron  conductivity  either  classical  or  anomalous  depending  on  plasma  parameters. 
Therefore,  the  self-consistent  electric  field  is  strongly  affected  by  density  gradients  and  anomalous 
transport. 

Initially,  one-dimensional  simulations  were  performed,  with  only  the  radial  (x)  direction 
numerically  resolved.  Without  azimuthal  direction  being  resolved,  2D  modes  like  Simon-Hoh 
cannot  be  excited  and  cannot  contribute  to  the  anomalous  transport  [Error!  Bookmark  not 
defined.].  The  Hall  parameter  values  were  about  an  order  of  magnitude  larger  than  observed  in 
the  experiment.  In  addition  to  scanning  the  gas  pressure,  parameters  sweeps  of  the  ion  mass  and 
magnetic  field  were  performed  and  allowed  us  to  identify  the  mode  responsible  for  the  anomalous 
transport  in  the  1-D  case  as  the  lower-hybrid  mode  [Error!  Bookmark  not  defined.]. 

In  2D  simulations  (see  Fig.7),  both  the  radial  and  azimuthal  components  are  resolved,  which 
allows  many  more  modes  to  get  excited  and  contribute  to  the  anomalous  transport.  Even  with  the 
reduced  device  size,  the  two-dimensional  simulations  are  somewhat  large.  500  x  500  cells  are  used 
to  resolve  the  Debye  sphere  in  the  plane  orthogonal  to  the  axial  magnetic  field  and  the  injected 
electron  beam,  and  a  time  step  of  5ps  is  used  to  resolve  the  plasma  oscillation.  25  macro  particles 
of  each  species  per  cell  were  found  to  be  sufficient  (quadrupling  the  number  did  not  significantly 
change  the  result).  With  these  parameters,  simulating  5us  (one  million  time  steps)  takes  just  over 
three  days  on  256  processor  cores  on  a  cluster  with  Intel  Xeon  processors  connected  with 
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InfiniBand.  To  date,  Penning-discharge  simulations  have  been  performed  using  over  370,000  core¬ 
hours,  equivalent  to  two  months  of  run  time  on  256  cores.  The  primary  parameters  scanned  have 
been  the  ion  mass,  gas  pressure  and  magnetic  field  [Error!  Bookmark  not  defined.].  The  physics 
is  more  complicated  in  two  dimensions  and  will  require  more  analysis,  but  the  mass  and  pressure 
scaling  of  the  Hall  parameter  seem  to  be  in  good  agreement  with  experiment  as  shown  in  Fig.  6. 

The  spoke  dynamics  is  complicated  and  it  is  difficult  to  determine  a  single  spoke  frequency. 
The  frequency  spectrum  is  very  broad  and  has  a  radial  dependence.  The  two-dimensional 
simulations  exhibit  rich  physics  that  need  more  theoretical  analysis.  Excitation  of  anti-drift  and 
low  hybrid  modes  together  makes  analysis  complex 


Figure  7.  Density  profile  (left)  is  peaked  in  the  center  and  is  different  for  different  ion  masses  (He,  Ar, 
Xe),  similarly  to  experimentally  observed.  Current  streamlines  on  top  of  potential  contours  (middle)  and 
electron-pressure  contours  (right)  at  2  ps. 

Experiments  on  Small  Low  Power  Magnetron  Discharges 

The  magnetron  discharge  studied  at  Stanford  is  shown  schematically  in  Fig.  8.  Shown  in  the 
figure  is  the  configuration  that  generates  a  5  mm  diameter  toroidal  discharge  described  in  several 
publications  [20  -  23].  As  part  of  this  program,  we  have  also  constructed  a  19  mm  diameter  plasma 


of  a  similar  construction  [24]  (see  photo  in  Fig. 
9).  In  this  configuration,  the  plasma  forms 
between  a  graphite  cathode  and  a  transparent 
indium  tin  oxide  anode  providing  direct  optical 
access  for  emission  and  laser-induced 
fluorescence  studies.  In  the  5  mm  diameter 
discharge  variant,  a  17  mm  diameter  Co-Sm 
ring  magnet  with  an  iron  core  generates  the 
magnetic  field  topology.  The  maximum  field  at 
the  cathode  is  ~1  T.  Axial  gradients  are 
generated  in  both  the  plasma  density  and 
magnetic  field.  These  gradients  drive  rotating 
structures  that  are  attributed  to  resistive  drift 
instabilities  (see  theory  below).  The  frequency 
that  we  see  in  these  discharges  is  much  higher 
than  the  rotating  structures  seen  in  the  Penning 
discharges  and  the  mode  structures  depend  on 
operating  conditions. 


Figure  8.  Schematic  of  5  mm  diameter 
magnetron  discharge  plasma.  Insets  show  high 
speed  camera  frames  of  emission  recorded 
through  the  ITO  anode,  and  side  photos  taken 
with  a  standard  DSLR  camera. 
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Figure  9.  DSLR  side-view  photograph  of  the  19 
mm  diameter  discharge.  The  left  is  the  case  of  a 
relatively  small  (~4  mm  gap),  whereas  the  right  is 
a  >  20  mm  gap.  The  discharge  shown  here  is 
operating  on  xenon. 


274  V,  7.6  mA 


We  have  developed  a  theory  for  the  rotating 
structures,  with  the  goal  of  predicting  the 
frequency  dependence  with  discharge  voltage 
(E-field)  and  the  mode  transitions  seen  with 
varying  discharge  voltage.  As  described  below 
the  theory  treats  the  disturbances  as  drift  waves 
driven  by  axial  gradients  in  both  density  and 
magnetic  field.  The  structures  are  well  defined 
and  very  coherent  (see  Fig.  10),  rotating  at 
speeds  of  2  -  4  x  103  m/s  -  considerably  lower 
than  the  local  drift  velocity  (~5  x  104  m/s).  The 
direction  of  these  coherent  modes  all  seem  to  be 
in  the  -  E  x  B  direction.  We  note  however  that 
the  direction  of  the  external  field  may  not 
necessarily  be  the  direction  of  the  local  field,  as 
we  discuss  below. 

The  mode  structure  seems  to  be  very 
sensitive  to  variations  in  the  discharge  voltage, 
undergoing  very  abrupt  transitions.  Within  a 
particular  mode,  the  mode  frequency  appears  to 
be  inversely  dependent  on  discharge  voltage 
(see  Fig.  11,  left  frame).  Increasing  discharge 
voltage  results  in  transitions  toward  higher 
mode  number  (shorter  wavelength)  as  can  be 
seen  in  Figs.  10  and  11.  Our  theoretical  studies 
suggest  that  the  mode  transitions  are  due  to  the 
strong  wavelength-dependence  of  the  cut-off  in 
the  growth  rate  of  this  gradient  drift  instability. 
Analysis  carried  out  with  segmented  ITO 
anodes  confirm  the  rotational  structure  seen  in 
the  high  speed  video  (see  Fig.  11,  right  frame). 
A  wavelet  dispersion  analysis  confirms  this 
strong  coherence,  with  a  resonant  frequency  at 
250  kHz  and  a  wavenumber  (1/a)  of  200  m"1. 
A  spectral  power  density  verses  frequency 
extracted  from  any  of  the  anode  segments 
shows  activity  between  10  and  50  MHz.  A  bi¬ 
coherence  analysis  suggests  the  possibility  of 
wave- wave  interactions  between  15  and  30  kHz 
[24] .  It  is  noteworthy  that  although  the  current 
in  each  segment  fluctuates  strongly,  as  shown  in  the  right  frame  of  Fig.  1 1,  the  sum  of  the  currents 
through  the  three  segments  as  well  as  the  remaining  anode  segment  (see  inset  in  Fig.  11 -right) 
results  in  a  nearly  non-fluctuating  total  current.  This  implies  that  these  structures  conduct  the  entire 
discharge  current.  The  discharge  current  is  anomalous  (i.e.,  cannot  be  accounted  for  by  classical 
scattering  collisions)  and  so,  as  was  found  in  the  studies  of  the  Penning  discharge,  implies  that 
these  structures  are  responsible  for  driving  anomalous  cross  field  transport.  The  predicted 
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Figure  10.  Frames  from  a  high  speed  camera  of 
emission  through  the  ITO  anode  depicting  rotating 
structures  and  their  mode-dependence  on  discharge 
voltage. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


frequencies  are  in  good  agreement  with  experimental  trends  (Fig.  11).  The  theory  (discussed 


Figure  11.  (left)  Measured  and  predicted  rotation  frequency  vs.  voltage,  (right)  Segmented  anode  signals 
during  a  strong  m  =  3  azimuthal  mode. 

below)  predicts  the  frequencies  to  scale  as  /  oc  (E0  /  mf12 ,  similar  to  that  suggested  by  the  studies 
of  Sakawa  et  al  [12]. 

Interesting  features  are  observed  when  the  voltage  is  set  such  that  the  system  resides  in  a 
metastable  state  between  modes  m  =  3  and  m  =  4,  or  between  m  =  4  and  m  =  5.  In  these  metastable 
states,  the  signal  in  any  single  anode  segment  fluctuates  in  a  somewhat  random  manner,  as  depicted 
in  the  left  frame  of  Fig.  12.  The  fluctuation  seems  to  be  quasi-turbulent.  High  speed  images  reveal 
considerable  absence  of  coherence,  and  the  segmented  anode,  which  serves  also  as  a  probe  of  the 
local  current  fluctuation,  show  a  rich  blend  of  frequencies,  with  fluctuations  at  very  high 
frequencies,  well  above  the  frequency  of  the  coherent  state.  The  panel  on  the  right  side  of  Fig.  12 
shows  a  spectral  power  distribution.  Note  strong  frequencies  as  high  as  3.5  MHz,  with  nearly  all 


Figure  12.  (left)  Segmented  anode  signals  when  state  is  between  m  =  4  and  m  =  5.  (right)  Corresponding 
power  spectral  plot  with  five  dominant  frequencies  indicated. 
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of  the  frequencies  higher  than  the  frequencies  of  the  coherent  modes  for  the  adjacent  m  =  4  or  m 
=  5  disturbances.  The  higher  frequencies  suggest  a  non-linear  mechanism  for  energy 
redistribution. 

A  simple  assignment  of  the  main  frequency  components  suggests  a  non-linear  wave-wave 
interaction  typical  of  a  turbulent  energy  cascade.  For  example,  for  the  results  shown,  we  see  that 
there  are  five  dominant  frequencies  associated  with  this  “turbulence”.  We  see  by  inspection  that 
to  a  great  extent, /1+/2  =fi,  andfi+f}  =fs,  suggesting  several  three-wave  coupling  possibilities.  The 
wave  dispersion  map  (see  Fig.  13)  from  the  current  that  is  collected  by  the  three  segments  also 

suggest  that  the  wavenumbers  of  these 
interactions  confirm  a  possible  three  wave 
coupling,  although  the  size  of  the  segments  (as 
seen  from  the  inset  of  the  left  frame  in  Fig.  11) 
preclude  this  assignment  with  certainty.  We 
also  investigated  the  characteristics  of  the 
discharge  over  higher  frequency  ranges  (up  to 
250  MHz)  and  found  between  10-40  MHz  a  fast 
mode  propagating  in  the  negative  ExB 
direction  with  a  phase  velocity  of  ~2xl 06  m/s. 
A  wavelet  analysis  indicated  the  existence  of 
high  turbulence  intensities  near  the  main  wave 
propagation  modes  indicative  of  energy 
dissipation  at  these  frequencies.  This  fast  wave 
seems  to  propagate  independently  of  the  lower 
frequency  modes,  whether  these  modes  are 
well  organized  or  not. 

Theory  of  Drift  Waves  in  Small  Low  Power  Magnetron  Discharges 

Because  of  the  relatively  high  pressure  (~  10-20  Pa)  and  plasma  densities  (~1013cm-3),  the 
plasma  in  these  magnetron  discharges  is  quasi-neutral.  The  millimeter  scale  produces  strong 
plasma  and  magnetic  field  gradients  (~1014  cm-4  and  5  T/cm,  respectively),  and,  as  a  consequence, 
we  interpret  our  results  using  models  for  gradient-driven  drift  instabilities  originally  developed  to 
explain  low  frequency  waves  seen  in  Hall  thrusters  [25,  26],  We  extend  the  model  of  Ref.  26  to 
include  transport  along  the  magnetic  field,  expanding  the  domain  of  instability.  The  model  results 
are  in  good  agreement  with  the  observed  experimental  behavior:  waves  propagate  in  the  -E  xB 
direction  with  frequencies  between  100  and  500  kHz  exhibiting  strong  dependency  on  the  local 
electric  field.  To  obtain  this  agreement,  we  must  assume  a  field  reversal  within  the  discharge, 
which  is  plausible  when  density  gradients  are  strong  and  diffusion  drives  more  anode-bound 
electrons  than  demanded  by  the  external  circuit.  Field  reversals  are  commonly  seen  in  hybrid 
simulations  of  Hall  thrusters  (in  the  anode  region),  and  one  was  recently  hypothesized  to  explain 
electron  heating  and  unsteady  behavior  [27]  and  enhanced  ion  back-flow  [28]  in  HiPIMS  sources. 

The  analysis  considers  a  two-dimensional,  two-fluid  model  where  the  ions  are  non-magnetized 
and  electrons  are  strongly  magnetized.  We  unfold  the  annular  geometry  onto  a  Cartesian 
coordinate  system  with  an  electric  field  E0  oriented  along  the  discharge  axis  (x-direction,  positive 
from  anode  to  cathode),  and  a  magnetic  field  B0  taken  primarily  along  the  z-direction  (radially 
outwards).  The  background  plasma  (assumed  uniform  in  the  y-direction)  contains  plasma  density 
(n0)  and  magnetic  field  gradients  along  the  x-direction,  as  illustrated  qualitatively  in  Fig.  14.  We 
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Figure  13.  Wave  dispersion  map  associated  with 
the  fluctuations  captured  by  segmented  anode  when 
the  discharge  state  is  between  coherent  modes  m  = 
4  and  m  =  5. 
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presume  that  the  strongest  ionization  and 
emission  occur  just  beyond  the  cathode  fall, 
where  plasma  density  gradients  near  the  anode 
produce  a  field  reversal,  E()  <  0,  driving  ions 
towards  the  anode.  An  analogous  reverse  ion 
migration  was  measured  by  laser  induced 
fluorescence  in  the  near  anode  region  of  Hall 
thrusters  [29],  highlighting  one  of  the  possible 
the  similarities  between  these  DC  magnetrons 
and  Hall  discharge  thrusters. 

Linearized  species  continuity  and 
momentum  equations  for  the  plasma 
constituent  number  densities  and  velocities  are 
subject  to  Fourier  perturbations  to  describe  the 
dispersion  conditions  for  wave  propagation.  The  quiescent  plasma  electrons  experience 
diamagnetic  (isothermal)  and  E  x  B  drifts,  respectively,  that  lead  to  the  following  equation  for  the 
plasma  density  [30]: 

8n'  +  ■  V«,  -  2«,  (v,0  +  v,£„  )•  V  In  B„  +  ^  =  0 


Figure  14.  Illustration  of  the  assumed  plasma 
structure  within  the  discharge  for  the  gradient-drift 
model. 


dt 


dt 


This  equation  differs  from  that  in  Ref.  26  by  the  last  term,  which  describes  diffusion  along  the 
magnetic  field  lines.  As  in  Ref.  30,  we  assume  electron/ion  collisionality  along  the  5-parallel  (z- 
direction)  only.  Invoking  quasi-neutrality,  the  electrons  and  ions  diffuse  together  with  a  z-directed 
flux,r  ,  =  -DAdnedz .  Here,  Da  is  the  ambipolar  diffusion  coefficient. 

By  carrying  out  the  usual  perturbation  analysis  and  linearizing  the  equations,  we  find  the 
following  dispersion  relation  for  the  frequency  (to)  dependence  on  wavenumber  components 
perpendicular  to  B0  (ki)  and  parallel  to  E0  (kx): 
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co  +  klv ?  + 
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=  0 
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co  -co D 


Here,  cs=  ( kTJtm)m  is  the  ion  acoustic  speed,  Vi0  is  the  ion  drift  velocity  due  to  the  external  electric 
field,  vc  =  k2zDA is  the  collisional  dissipation  rate,  co  =  kyveD  ,  coa  =  kyveExB  and coD  =  -2kT,  / eBaLB  . 

In  the  absence  of  the  diffusive  term,  the  dispersion  will  have  identical  characteristics  to  that  derived 
in  Ref.  26,  and,  with  E0  (and  Vi0)  <  0,  the  region  of  instability  has  a  distinct  long  wavelength  cutoff 
that  depends  on  background  plasma  conditions.  We  calculate  the  predicted  growth  rate  of  the 
unstable  branch  (imaginary  component  of  the  dispersion  solution  assuming  real  wavenumbers) 
using  experimental  conditions,  E0=~l  1.6  kV/m,  Ln=2  mm  (E0Ln=-23.l  V),  Lb=0.5  mm,  kTc 
=3.1  eV,  Bo=0.6  T,  kx  =7.0  mm-1,  kz= 2.5  mm-1,  and  Vi0  =  -933  m/s,  and  plot  the  result  as  the  solid 
line  in  Fig.  15(top).  The  maximum  growth  rate  is  somewhere  between  m  =  3  and  m  =  4,  and 
increasing  the  electric  field  shifts  the  cutoff  and  peak  towards  longer  wavelengths.  Including  the 
diffusive  term  (using  Vi  =  5  *  107  s-1)  modifies  the  growth  rate,  extending  the  region  of  instability 
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beyond  the  original  cutoff  (see  the  dotted  line 
in  Fig.  15(a)).  As  the  discharge  voltage 
increases,  we  expect  an  active  m  =  4  mode  until 
the  peak  growth  rate  passes  a  critical  value, 
triggering  a  mode  transition  to  m  =  3.  The 
variation  in  the  peak  growth  rate  with  electric 
field  is  shown  in  Fig.  15(b).  When  increasing 
the  field,  we  assume  a  simple  scaling  that 
kT£  ~  E0  and  Vi0~£'01/2. 

The  real  component  of  the  frequencies  of 
the  unstable  root  tends  to  decrease  with 
increased  electric  field  while  within  the  range 
of  instability  for  any  single  mode.  As  illustrated 
by  the  solid  lines  in  Fig.  1 1  (left  frame),  once  a 
new  mode  is  preferred  (higher  m  for  voltage 
increases),  a  jump  in  frequency  is  predicted.  To 
compare  the  experiments  with  the  model,  we 
have  assumed  that  a  12  V  voltage  “hump” 
(E0Ln)  is  generated  when  the  discharge  voltage 
is  261  V.  No  other  adjustments  are  made  to  the 
predicted  dispersion.  The  values  chosen  for 
various  plasma  parameters  are  the  same  as 
those  listed  above  and  used  in  Fig.  15. 
Agreement  with  experiments  is  good, 
considering  the  uncertainties  in  the  properties 
used  as  base  conditions.  A  model  tested  without 
including  the  diffusive  terms  fails  to  achieve 
such  an  agreement. 

The  low  power  DC  magnetron  studied  as 
part  of  this  project  is  an  excellent  test  bed  for  understanding  ordered  (coherent)  structures,  which 
we  believe  can  be  attributed  to  gradient-driven  instabilities  in  cross-field  discharges.  The  observed 
mode  transitions  are  distinct,  coherent,  and  reproducible.  A  simple  theory  for  the  dispersion  of 
these  gradient-driven  drift  waves  seems  to  describe  these  instabilities  quite  well.  The  success  of 
the  model  rests  heavily  on  the  presumption  of  a  field  reversal,  generating  a  voltage  “hump”  that 
drives  a  reverse  (anode-directed)  ion  flow,  and  future  experiments  will  seek  confirmation  of  such 
a  structure.  Our  anticipation  is  that  advanced  diagnostics  will  enable  validating  this  field  reversal, 
and  hence  the  essential  physics  needed  to  drive  the  rotation  of  these  structures  in  a  direction 
opposite  to  the  imposed  electric  field. 

One  note  in  point  is  that  the  frequencies  of  these  structures  are  nearly  -1MHz.  We  believe  that 
these  structures  can  be  driven  to  higher  frequencies  with  external  forcing.  Their  external  control 
can  enable  new  applications  such  as  in  the  control  of  electromagnetic  waves  in  photonic  crystal 
waveguide  filters. 


m  =  kR. 
y  o 


Figure  15.  (a)  Growth  rate  of  the  instability  versus 
mode  number  for  the  diffusive  (dashed)  and  non- 
diffusive  (solid)  cases,  (b)  Growth  rate  variation  with 
mode  number  computed  for  a  range  of  values  of  the 
voltage  “hump”. 
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